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Abstract. The HST/ACS Survey of Galactic globular clusters (GGCs) is a HST Treasury 
project aimed at obtaining high precision photometry in a large sample of globular clusters. 
The homogeneous photometric catalogs that has been obtained from these data by Anderson 
et al. (2008) represents a golden mine for a lot of astrophysical studies. 
In this paper we used the catalog to analyse the properties of MS-MS binary systems from 
a sample of fifty GGCs. We measured the fraction of binaries (divided in different groups), 
studied their radial distribution and constrained the mass ratio distribution. We investigated 
possible relations between the fraction of binaries and the main parameters of their host 
GGCs. 

We found a significant anti-correlation between the binary fraction in a cluster and its abso- 
lute luminosity (mass). 
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1. Introduction 

Knowledge of the binary frequency in Globular 
Clusters (GCs) is of foundamental importance 
for a lot of astrophisical studies. 

Binaries play an important role in the dy- 
namical evolution of a clusters. Interactions 
with hard binaries pump kinetic energy into 
the cluster core, slowing the core collaps and, 
eventually, causing the core to reexpand, if the 
number of binaries is large enough. In general, 
binaries are a foundamental ingredient in any 
dynamical evolution model of a GC. 

Exotic stellar objects, like Blue Stragglers, 
cataclismic variables, millisecond pulsars and 
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low mass X ray binaries are believed to rep- 
resent evolutionary stages of close binary sys- 
tem. The determination of the fraction of bi- 
naries plays a foundamental step towards the 
understanding of the evolution of these pecu- 
liar objects. Furthermore, binary stars intro- 
duce systematic errors in the determination of 
the main sequence (MS) fiducial line and move 
it toward red colors with respect to its coiTect 
position. 

Finally, a correct determination of the mass 
and luminosity functions requires a correct 
measure of the fraction of binaries. 

Up to now, three main techniques have 
been used to measure the fraction ob binaries 
in GGCs (Hut et al. 1992). 
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The first one identifies binaries by measur- 
ing their radial velocity variation (eg. Latham 
1996). This method relies with the detection of 
each individual binary system but, due to the 
limits in sensibility of spectroscopy, these stud- 
ies are possible only for the brightest GGCs 
stars. 

The second tecnique is based on the search 
for photometric variables (eg. Mateo 1996). As 
well as the previous one, it is able to infer spe- 
cific properties of each binary system (like the 
measure of orbital period, mass ratio, orbital 
inclination). Unfortunatly, it is biased towards 
binaries with short periods and large orbital in- 
cUnation. Moreover these tecniques have a low 
discovery efficency and are very expensive in 
terms of telescope time because it is necessary 
to repeat measures in time. 

A thirth approach, that is based on the anal- 
ysis of the number of stars located on the red 
side of the main sequence (MS) ridge hne 
(MSRL) may represent a more efficient method 
to measure the fraction of binaries in a cluster 
for several reasons: 

- availabihty of a large number (thousands) 
of stars makes it a statistically roboust 
method; 

- it is cheep in terms of observational time: 
two filters are enough for detecting binaries 
and repeated measuremets are not needed. 

- it is sensitive to binaries with any orbital 
period and inclination 

This approach have been used by other groups 
( see SoUima et al. 2007 and references within) 
to study the population of binaries in GGCs. 
The relative small number of clusters that have 
been anahzed is consequence of the intrinsic 
difficulties of the method: 

- high photometric quality is required; 

- in some cases, the differential reddening 
spreads the MS and makes it more difficult 
to isolate the binary sequence; 

- an accurate analysis of photometric er- 
rors as well as a correct estimate of field 
contamination are necessary to disantangle 
real binaries from bad photometry and field 
stars. 



In this paper, we analyse the catalogs ob- 
tained by Anderson et al. (2008) from HST 
ACS/WFC data. We exploited both the home- 
geneity of this dataset, and the high photomet- 
ric accuracy of the measures to derive the frac- 
tion of binaries in the central regions of fifty 
GGCs. 

2. Outliers in the Color Magnitude 
Diagram 

Binaries that are able to survive in the dense 
enviroment of a globular cluster are so close 
that even the Hubble Space Telescope (HST) 
is not able to separate their single components. 
For this reason, light coming from each star 
will combine, and the binary system will ap- 
pear as a single pointlike source. 

In this paper we will take advantage from 
this instrumental limit to search for binaries by 
analising their peculiar position in the CMD. 

In the most general case, if we consider 
two stars in a binary system and indicate with 
nil, ni2, F\ and F2 their respective magnitudes 
and fluxes, a simple algebric count demon- 
strate that the binary will appear as a single ob- 
ject with a magnitude: 

mbin ^ nil -2.5 log(l + jr) 

In the case of a binary formed by two MS stars 
(MS-MS), fluxes are related to stellar masses 
(Ml, M2), and its luminosity depends on the 
mass ratio q = M2/M1 (in the following we 
will assume M2 < Mi, q < 1). The bina- 
ries formed by an equal mass pair form a se- 
quence parallel to MS, and ~ 0.75 magnitudes 
brighter. When the masses of the two compo- 
nents are different, the binary will appear red- 
der and brighter than the MS and it will be lo- 
cated in a CMD region on the red side of the 
MSRL. 

An obious consequence of this analysis is 
that our capabiUty in detecting binaries de- 
pends meanly by the photometric quality of 
the data. Binaries with large mass ratios have 
a large distance from MSRL and are relatively 
easy to be detected. On the contrary, a small 
mass ratio pushes binaries near the MSRL and 
makes it hard to separate them from single MS 
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stars. Moreover, the poorer photometry of faint 
stars hmits the luminosity (mass) range where 
they can be detected and studied. 

3. Method 

The Hmited photometric precision makes it im- 
possible to measure the overall population of 
binaries even in a small region of clusters. For 
this reason, in this paper, we do not pretend 
to measure the global fraction of binaries in 
a cluster, but will limit our study to particu- 
lar subsamples of them. Each group is formed 
by objects that share all the same properties in 
terms of luminosity and mass ratio. 

We isolated three samples of high mass ra- 
tio binaries (defined as the binary systems with, 
q > 0.5, 0.6 and 0.7) and separately studied the 
properties of each group. In addition we de- 
rived also the global fraction of binaries. 

We performed our study in the magnitude 
range 3.75 < A{Ifsi4w) < 0.75 below the main 
sequence turn off. The extremes of this interval 
will be indicated with Ibright and I/aint- 

3.1. High q binary fraction 

In order to measure the fraction of high q bina- 
ries, we divided the CMD in two regions (see 
FiglB: 

A region (A) that includes all the MS single 
stars and the binaries with a primary star with 
hright > I > I faint- It is formed by three subre- 
gions. The first one (Ai) includes all the MS 
single stars and MS-MS binaries with small 
mass ratios; it is limited by dashed lines in 
Fig[T]and corresponds to the CMD portion with 
a color distance from the MS ridge line smaller 
than three times the MS dispersion; the second 
(A2) includes all the binary candidates with 
high mass ratios. In Fig[T]it corresponds to the 
CMD portion on the red side of the A i region 
and is delimited by; 

- the track formed by a binary system with 
a primary star with / = hright and a mass 
ratio ranging from to 1 on the top; 

- the corresponding track for a binary system 
with a primary star of mass //«,„, on the bot- 
tom; 
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Fig. 1. Grey areas are the regions of the 
NGC 6121 CMD adopted to select all the (sin- 
gle and binary) cluster stars (right) and the can- 
didate binaries with q>0.6 (left) in a range of 3 
h%uw magnitudes. 



- the ridge line for an equal mass binary sys- 
tem on the red side. 

The third region (A3) contains all the binaries 
with q ~ \ that are shifted by photometric er- 
rors to the right of an equal mass binaries fidu- 
cial line. It is adiajent to the region A2 and it is 
limitated by the ridge line for an equal mass bi- 
nary system shifted in color by three times the 
main sequence photometric dispersion on the 
left side. 

The second region B is defined as the por- 
tion of the region (A) on the red side of the 
track formed by a binary star with q - q,,- and 
it includes all the binaries formed by a primary 
star with Itright > I > I/aint and a mass ra- 
tios greater than a threashold value (qtr)- In this 
work we separately studied the samples of bi- 
naries with mass ratios larger than ^,^=0.5, 0.6 
and 0.7. 

Unfortunatly, regions A and B are popu- 
lated by field stars, while chance superposition 
of two unrelated stars (apperent binaries) may 
reproduce the behaviour of a genuine binary 
system. 
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To estimate the quantity of back- 
ground/foreground objects that casually 
overlap the cluster CMD we used the galactic 
models of Girardi et al. (2008) (with the 
ecception of seven clusters, where we could 
isolate field stars through proper motions). 
The fraction of apperent binaries has been 
quantified by performing artificial star tests. 

We further applied a the tecnique described 
in Milone et al. 2008 to correct the spread in 
color caused by differential reddening and/or 
spatially dependent zero point photometric er- 
rors. 

In order to measure the fraction of high 
mass ratio q binaries, we started by deriving the 
observed numbers of stars in regions A {Ngg^) 
and B {Ng^^). Then we evaluated the corre- 
sponding values of arificial stars {N^jfj and 
A^^^^) and field stars (Njj^^jj and A^^^^^^,) The 
correct numbers of real, field and artificial stars 
are calculated as N - Yj where c,- is the 
completeness. 

High mass binary fraction is calculated as 

f _ ^^OBS FIELD _ ^^ART 

J bin — mA _fjA fjA 

'^OSS '^FIELD ^^ART 

3.2. The global binary fraction 

In the case of MS-MS binaries with small mass 
ratios, the distance from the MSRL is compa- 
rable to the size of photometric errors in color. 
Therefore these objects appear mixed to single 
stars, and a more sophisticate statistical anal- 
ysis is required to derive their contribution to 
the global fraction of binaries. 

In order to estimate the global fraction of 
binary systems we adopted a statistical method 
which is based on the comparison of the ob- 
served data with more then 10,000 simulated 
CMDs enriched by different fractions of bina- 
ries with a given fiq). Details of this tecnique 
are outside the purpose of this work and will 
be expleined in a farcoming paper (Milone et 
al. 2008). We want emphasize here that the 
adopted statistical approach is very sensitive 
to small inaccuracies and the results that we 
present are to be considered just as approxi- 
mate estimates of the real fraction of binaries 



4. Results 

We investigated if the fraction of binaries de- 
pends on the radial distance from the cluster 
center. To this aim we divided the ACS field 
into four concentric annuli, each containing 
roughly the same number of stars, and mea- 
sured the fraction of binaries in each of them. 

In (at least) the ~ 60 % of the 37 GGCs 
with good photometry in the innermost re- 
gions, binaries are more centrally concentrated 
than single MS stars. In the other objects, their 
distribution is consistent with a flat distribu- 
tion, probably, as a consequence of the small 
radial coverage of our field of view. 

In the following, we present preliminary 
results that involve binaries selected from dif- 
ferent cluster regions: 

- inside the core radius (Icore sample); 

- between the core and the half-mass radius 
(HM sample); 

- outside the half-mass radius {Ohm sample). 

Even if data used in this paper are homoge- 
neous as they came from the same observ- 
ing facility (ACS/HST) and have been reduced 
adopting the same tecniques, their photometric 
quality (and completeness) may vary from one 
cluster to the other, mainly, as a consequence 
of the different stellar densities. 

For this reason, it was possible to include 
in the Icore sample only 35 out 50 GGCs. In 
addition, the limited ACS field of view reduced 
the number of GGCs with HM and Ohm sam- 
ples to 46 and 29 respectively. 

We explored possible relations between the 
fraction of binaries and the main parameters 
of their host GGCs absolute visual magnitude, 
metallicity, collisional parameter, core and half 
mass relaxation time, central density and con- 
centration (from Harris et al. 1996, 2003). 

We found an highly significant anticorre- 
lation between the binary fraction and the total 
cluster luminosity with clusters with fainter ab- 
solute luminosity have higher binary fractions. 
This anti-correlation is shown in Fig|2]for the 
Icore sample and in Figs[3] and for the HM 
(bottom panel) and Ohm (upper panel) sam- 
ples. 

A similar anticorrelation has been found by 
Piotto et al (2004) and Moretti et al. (2008) 
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Fig. 2. Fraction of binaries with q > 0.5, 0.6 

and 0.7 and global fraction of binaries in the 
core as a function of the host GGCs luminosity. 
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Fig. 3. Fraction of binaries with q > 0.5 for 
the Ohm (top) and HM samples (bottom) 



between the frequency of blue stragglers and 
the luminosity of the cluster. It is very tempt- 
ing to connect BSS and binaries populations 
in GGCs: the BSS frequency can be in fact 
related to the evolution of the binary fraction 
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Fig. 4. Fraction of binaries with ^ > 0.5 in the 
core as a function of the collisional parameter 
(Davies et al. 2004) 



due to encounters as described by Davies et al. 
(2004). We found only a marginal correlation 
between the binary fraction and the cluster col- 
lisional parameter 
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